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Abstract. - Using molecular dynamics simulations we investigate the relaxation dynamics of a 
supercooled liquid close to a rough as well as close to a smooth wall. For the former situation the 
relaxation times increase strongly with decreasing distance from the wall whereas in the second 
case they strongly decrease. We use this dependence to extract various dynamical length scales 
and show that they grow with decreasing temperature. By calculating the frequency dependent 
average susceptibility of such confined systems we show that the experimental interpretation 
of such data is very difficult. 



Motivation. - The details of the mechanism giving rise to the dramatic slowing down 
of the dynamics of glass- forming liquids upon supercooling are still unknown (see, e.g. 
Although the mode-coupling theory of the glass transition allows to rationalize many features 
of the relaxation dynamics of these systems 0] , the answers to certain important questions 
(e.g. the relaxation dynamics at low temperatures) are still unknown. A further popular 
approach is the phenomenological concept of "cooperativity" , introduced by Kauzmann fl , 
and Adam and Gibbs Q . A typical example of this cooperativity is the so-called "cage-effect" , 
i.e. the fact that in a dense liquid each particle is surrounded by neighboring particles which 
form a temporary cage. In order to allow the particle to change its position the cage has to 
open up. However, each of the particles of the cage is itself also caged and hence can move 
only if other particles make room. Therefore one can conclude that the particle motion is 
collective and there exist "cooperatively rearranging regions" (CRR's) within the liquid. The 
typical size of a CRR is postulated to grow with decreasing temperature, hence "rationalizing" 
the slowing down of the dynamics 

Experimentally it is difficult to test the concept of the CRR's since usually one does not 
have direct access to the dynamics of single particles. Therefore many studies have focused 
on investigating systems in spatial confinement. If CRR's do exist and grow with decreasing 
temperature the dynamics should differ from the bulk behavior as soon as the size of the 
CRR's at a given temperature becomes comparable to the system size. Indeed, almost all 
experiments on glass formers confined to porous host material J6j-|l2| and supported (or even 
free standing) films [p~3| |lT|| do indeed show a relaxation dynamics that differs from the one 
in the bulk. However, so far it has not been possible to give a conclusive interpretation of 
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experimental results, since, e.g., one sometimes finds that the dynamics in confined systems 
is faster than the one of the bulk, whereas in other systems it is slower (see, e.g., |]s|, |l(|) . 
The main reason for the diversity in experimental findings (see, e.g., the review article ]18[) 
is the influence of secondary effects beyond the picture of CRR's: On the one hand density 
effects can be expected to play a crucial role. Since in experiments only the average density 
is accessible, one can imagine a situation where the confined liquid shows strong local density 
variations, e.g. density oscillations due to layering effects, which will strongly influence the 
dynamics as well. Also the interaction between the surface and the liquid is most likely very 
important [ [T4] , |l5] |. If the liquid particles are sticking to the wall, their mobility is strongly 
suppressed and because of cooperativity also the dynamics of particles a certain distance away 
from the wall will be slowed down, which in turn will influence the dynamic properties of the 
whole sample. The opposite should hold in a situation where the liquid particles can slip along 
the boundary. 

The goal of the present paper is to use molecular dynamics computer simulations to inves- 
tigate how these boundary effects influence the relaxation dynamics of confined glass-forming 
liquids and to what extend growing length scales can be extracted in such systems. 

Simulation. The model liquid under investigation is a binary mixture of particles 
interacting via a Lennard- Jones (LJ) potential V a p{r) = ie a /3[(a a fi/r) 12 — (cr a fi/r) 6 ] with 
a, [3 £ {A,B}, cut-off radii 7^^=2.5 • <j a p and interaction parameter eaa = 1-0, £>aa = 1-0, 
£ab = 1-5, cab = 0.8, €bb = 0.5, and ctbb = 0.88. In the following we will use oaa and eAA 
as units of length and energy, respectively, setting Boltzmann's constant fcs = 1, and measure 
time in units of \J mcr AA /48eAA, where m is the mass of the particles. Previous simulation in 
the bulk have shown that at low T the dynamics of this mixture quickly slows down and that 
its mode-coupling temperature T c is around 0.435 C^). 

We considered two types of walls: a rough and a smooth one. The rough wall was realized 
by freezing a slice of thickness 2.5oaa of the LJ liquid and applying the same LJ interaction 
plus an additional hard core potential to prevent the liquid particles from penetrating into 
the wall. More details on this type of wall can be found in |£ll 

To mimic a smooth wall at a location zw we applied an external potential of the form 
V aW (z) = (4/45)7T / 9 W CT AB e QW (o'AB/(z - z w )) 9 (with e A w = 1-0, £bw = 3.0). By choosing 
walls at zw = —0.65 and zw = 15.65 we make sure that a film width of D = 15.0 is realized. 

The data for the rough wall presented here comes from a simulation of a film with area 
LxL and thickness D (L — 12.88, D = 15.0), where periodic boundary conditions are applied 
in the film plane. It contains 2400 A and 600 B particles giving an average density of 1.2 as 
in the bulk simulations in Ref . |lj]] . For the smooth wall the area of the film was four times 
larger. In order to improve the statistics we averaged the results over 16 independent samples. 
The equations of motion were integrated with the velocity form of the Verlet algorithm, using 
at low T a time step of 0.02, and the starting configurations for the microcanonical production 
runs were carefully equilibrated. 

As already mentioned in the previous section it is most important to avoid that the confin- 
ing walls change the structural properties of the enclosed fluid. This request is by construction 
fulfilled trivially in the case of the rough wall. For the smooth wall, however, we found that 
close to the wall strong layering effects occur [p2[ . To avoid this problem we modified the 
potential energy of the system by adding a term that coupled directly to the deviation from 
a constant density profile, i.e. configurations whose profile is not constant (within a certain 
fluctuation) are energetically disfavored. This modification does indeed allow to obtain a 
system with an essentially constant density profile and we have checked that its structure is 
indeed the same as the bulk one [^2| ) . 
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Fig. 1 - Structural relaxation time r q (z) at q = 7.2 as a function of particle distance z from the wall 
for (a) rough and (b) smooth surface at different temperatures. The large diamonds are the bulk 
values and the long dashed lines and the solid ones are fits according to Eqs.(m) and (0), respectively. 



Results. - We now discuss the relaxation dynamics of the system as a function of tem- 
perature and distance from the wall. Although in the following we focus on the self part of 
the intermediate scattering function, all the conclusions are also valid for other observables, 
such as mean squared displacements or the van Hove correlation functions 1 2^| . 

From earlier investigations of systems confined between rough surfaces [20, 21j we know that 
the local particle dynamics depends on the distance z from the wall. Therefore we introduce 
a generalization of the incoherent intermediate scattering function by f s (q, z,t) = N^ 1 X)j=i 
(exp [iq ■ (ij(t) — rj(0))] 5(zj(0) — z)), i.e. F B (q, z, t) considers only particles that at t = had 
a distance z from the wall pl[ . We will only present data for A particles and wave vectors 
q parallel to the wall with |q| = 7.2, the location of the maximum in the structure factor. 
Data for other values of |q| or the B particles looks qualitatively similar. At low temperatures 
F s (q, z, t) shows a two step relaxation. Hence we can characterize the a— relaxation time T q (z) 
byF B (q, 2 , [ 7v)= e - 1 @,||]. 

In Fig. [I] we show the relaxation times T q (z) as a function of z at different temperatures. 
For particles far away from the wall F s (q,z,t) shows bulk behavior and therefore the char- 
acteristic relaxation times coincide with the bulk values (filled diamonds). Approaching the 
rough surface the dynamics is slowed down dramatically, and r q (z) grows continuously over 
several decades in time. In contrast to this the dynamics of the system with smooth surfaces 
accelerates with decreasing z. Furthermore one can see that, in both situations, the region 
affected by the wall expands with decreasing temperature in that for high T the bulk behavior 
is realized already at z > 2.0, while at the lowest T the influence of the wall extends almost 
to the center of the film. In the following we will use this increase to define a length scale. 

In Ref. it was shown that close to the surface the z dependence of T q (z) can be described 
well by the empirical Ansatz 



r 9 (z) = /,(T) e x P [±A 9 (T)/(z-z p )] 1 



(1) 



with three free parameters f q (T), z p = 



—0.5 ±0.15 and A q (T), the latter weakly temperature 
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Fig. 2 - Different dynamic length scales (see text for definition) in confined systems as function of 
inverse temperature in a logarithmic plot. Comparison with static and dynamic length scales in bulk 
systems. Data for £ c i is taken from Ref. [0. 



dependent. (Here and in the following the positive (negative) sign corresponds to systems with 
rough (smooth) surfaces.) The z— range for which this fit works well increases with decreasing 
T (dashed curves in Fig. [j]). The location of the crossover to bulk behavior in the center (see 
data and fits for T = 0.6 in Fig. |l|a) can be used to define a characteristic length scale S whose 
T— dependence will be discussed below. 

An alternative Ansatz for the z— dependence of T q (z) is a function depending on exp(— z/£q), 
with a characteristic length scale £o(T). We find that the functional form 



In 



(r q (z)/T qtOC ) 



±i 



A(T) ■ exp [-z/to(T)] 



(2) 



does indeed describe the z— dependence for all values of z (solid lines in Fig. [T]). Here r 9j0 o is 
the relaxation time of the system in the bulk. 

The length scales z(T) and £o(T) are obtained directly from r q (z). In Ref. pit] we have 
shown that for the case of the rough surface the whole time and z— dependence of F s (q, z,t) 
can be described very well by the Ansatz 



F 8 (q, z, t) = F s bulk (q, t) ± a(t) exp - (*/£(*)) 



(3) 



The time dependence of the length £(t) is smooth and shows a maximum, thus showing that 
the influence of the wall on the dynamics is maximal on the time scale where this maximum 
occurs (which is on the order of the a— relaxation time of the bulk). The value of the maximum 
in thus allows to define a dynamical length scale £ ma x- 

In the case of smooth surfaces the Ansatz (||) is not very useful to define a length scale. 
Although the data is still described reasonably well by Eq. (||), £(t) is growing monotonically 
and therefore it is not possible to read off a £ max p2f . 
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Fig. g is an Arrhenius plot of the T— dependence of z, £o and £max and from this graph we 
see that these length scales grow like oc exp[E/T]. Note that the activation energy E « 1.1 
depends neither on the definition of the length scale nor the type of wall, which gives evidence 
that the length scales do indeed characterize a length scale intrinsic to the system. Within the 
T— range investigated the growth of this scale is rather small, e.g. £o grows only by a factor 
of 3.5 between T — 2.0 and T — 0.5 while bulk relaxation times increase by several orders of 
magnitude. Furthermore we find no evidence for a divergence at a finite temperature close to 
the investigated temperature region, in agreement with Ref. pl[ . 

Finally we compare these length scales with the ones that can be identified already in the 
bulk. Using the van Hove autocorrelation function one can identify the most mobile particles 
in the system (23) . It turns out that these particles are not distributed randomly but instead 
form temporary clusters (dynamical heterogenieties), the size of which grows with decreasing 
temperature (2^] . This size is included in Fig ^ as well and we find it to be comparable to the 
dynamic length scale from the present simulation if T > 0.5. 

As an example for a static length scale (in bulk and film systems) we consider the decay 
length of gAA{r), the radial distribution function for AA correlations. The envelope of <7aa(?") — 
1 is described well by exp(— r/£ s ), which thus defines a length scale £ s . This length scale shows 
only a very weak T-dependence (Fig.||), in contrast to the dynamical length scales. 

Interpretation and Comparison with Experiment. - Having shown how the nature of the 
walls influences the relaxation dynamics of the particles of the confined system, we now discuss 
the signature of the observed slowing down/ acceleration of the dynamics in an experiment in 
which the average dynamics of the particles is measured. This is motivated by the fact that in 
a real experiment it is usually not possible to determine the relaxation dynamics as a function 
of the distance from the wall whereas the average dynamics is directly accessible. 

Thus to compare our simulations with typical experimental data we have to calculate 
dynamic properties averaged over the whole system, such as F 8 (q, t), the integral of F„(q, z, t) 
over z. Although the single curves for F s (q, z, t) at different z look very similar for rough and 
smooth surfaces (2(], ^] and the z— dependence of the relaxation times r q (z) has the same 
functional form we find a qualitative difference in the averaged curves for i^ s (q, t). 

For systems with smooth surfaces the a-relaxation is described well by a stretched expo- 
nential law with a stretching exponent that is slightly lower than the corresponding bulk value, 
i.e. due to the superposition of different relaxation processes the curves are more stretched. 

In contrast to this, F s (q, t) for the rough walls shows a long time tail because of the 
huge relaxation times for particles at the surface [^0|. It is possible to describe the whole 
a— relaxation, i.e. also the mentioned tail, by the sum of two stretched exponentials, where 
the time scale of the slow "process" is about two orders of magnitude larger than the first 
"process" and the stretching is much more pronounced |2^| . Note however, that this is a 
purely phenomenological description of the data without any underlying physical motivation. 

Typically experiments on confined liquids probe the dynamics of the system by measuring 
the frequency dependence of various susceptibilities (light-and neutron scattering scattering 
experiments, dielectric measurements). Hence we have calculated the dynamic susceptibility 
x"(q, w) = w/(2fcBT)5 s (q, w), where S B (q,u) is the time-Fourier transform of F s (q, t). 

In Fig. |||a/b the frequency dependence of x"(q, u>) at a low T is shown for bulk systems 
as well as for films with smooth and rough surfaces. (Note that we show only the frequency 
range of the a— relaxation. The microscopic peak is around u = 1). Fig. |^a shows that a 
smooth surface has a broader a— peak than the bulk, in agreement with the observation that 
in the time domain the stretching exponent is smaller. Furthermore the position of the peak 
is shifted to slightly higher frequencies, which can be understood from the fact that close 
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Fig. 3 - a/b: Dynamic susceptibility Xs( c li tJ ) f° r a LJ liquid in the bulk and confined to a film 
with smooth and rough surfaces at q = 7.2, T = 0.5. The dotted lines are stretched exponentials, 
c/d: Experimental data for the imaginary part or the dielectric susceptibility of salol in the bulk and 
confined in pores. The dotted lines are fits with Havriliak-Negami functions. Adapted from Ref. fl9|. 



to the surface the particles move faster than in the bulk. In contrast to these rather small 
differences between the curve for the bulk and the one for the smooth surface, the spectrum 
for the system with the rough wall differs strongly from the one of the bulk (Fig. |^b). First 
of all we note that the location of the peak is shifted to smaller frequencies. More important 
is, however, that the curve for the rough surface shows a pronounced shoulder left to the a— 
peak. From Fig. |l| we know that it would be wrong to ascribe this shoulder to the presence 
of a second relaxation process, since the relaxation times are a smooth function of z. Instead 
this shoulder is just due to the superposition of a continuum of relaxation processes with very 
different relaxation times. Also included in the figure are the Fourier-transforms of the two 
stretched exponentials (dotted lines) which, noted above, describe the whole a— relaxation in 
the time domain. 

The same qualitative behavior of relaxation spectra for liquids in confinement is also seen 
in many experiments. As an example we take dielectric data of the simple glass former salol 
confined in Vycor glass from Ref. Q. In this "quasi" -van der Waals liquid H-bonds are mainly 
of intramolecular nature and therefore the interaction is to a first approximation comparable to 
the van der Waals system studied here. Samples with uncoated pores, where H-bonds between 
the molecules and the wall can form, correspond to the situation of a rough surface. If the pore 
surface is coated, i.e. the formation of the mentioned H-bonds is prevented, the interaction 
between the wall and the liquid becomes weak, i.e. a smooth surface is realized. In Fig. [|c/d 
we show the data from Ref. || for the bulk, as well as porous systems with a native surface 
and a treated one. We see that from a qualitative point of view these three spectra are very 
similar to the ones obtained in our simulation for the three different situations (broadening 
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of the peak in the confined system, presence of a shoulder at low uj in the untreated surface, 
etc.). (Similar experimental results can be found in Refs. Also note that the increase of 

the spectra at low u> is related to the Maxwell- Wagner polarization of the sample Q which has 
nothing to do with the structural relaxation of the system.) We emphasize, however, that for 
our system the presence of the second peak in the susceptibility is not related to the existence 
of a layer at the surface that relaxes orders of magnitudes slower than the rest of the liquid, 
a popular interpretation of such a feature |7[-p|, |l6|, fL7[ p4|| . Instead our analysis of the local 
dynamics has shown that this secondary peak is just the result of averaging over particles that 
have a continuous and monotonous distribution of relaxation times. 
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